Sub-Poissonian atom number fluctuations by three-body loss in mesoscopic ensembles 
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We show that three-body loss of trapped atoms leads to sub-Poissonian atom number fluctuations. 
We prepare hundreds of dense ultracold ensembles in an array of magnetic microtraps which undergo 
rapid three-body decay. The shot-to-shot fluctuations of the number of atoms per trap are sub- 
Poissonian, for ensembles comprising 50-300 atoms. The measured relative variance or Fano factor 
F — 0.53 ± 0.22 agrees very well with the prediction by an analytic theory (F = 3/5) and numerical 
calculations. These results will facilitate studies of quantum information science with mesoscopic 
ensembles. 

PACS numbers: 03.75.Be, 05.40.-a, 42.50.Lc, 67.85.-d 



The study and control of particle number fluctuations 
in ultracold atomic systems has revealed a rich variety of 
intriguing quantum phenomena [3,0)S0)B@|> an d offers 
the potential to boost performance in cold atom tech- 
nologies. Motivated largely by the prospects for quan- 
tum metrology @, B @] , several recent experiments have 
demonstrated the suppression of relative fluctuations be- 
tween small atomic samples distributed over two or more 
traps or internal states, leading to number difference or 
spin squeezing and entanglement 
contrast, however, work on suppressin 
ber fluctuations has been limited 14|, []__ 
cial, for example, in quantum information science using 



11, h 



130- By 

absolute num- 
. This is cru- 



mesoscopic atomic ensembles 16, 17, 3, HI, where re- 
cently observed collective excitations produced via Ry- 
dberg dipole blockade [2^, Hi], [22| could be exploited. 
Trapped ensembles would benefit from a V^V collective 
enhancement of the Rabi frequency over single atoms, al- 
lowing fast quantum operations. However, intrinsic atom 
number fluctuations would adversely affect the fidelity. 

In this Letter we show explicitly that three-body loss 
naturally reduces the shot-to-shot fluctuations of the ab- 
solute atom number in a trap to sub-Poissonian levels. In 
experiments, random particle loss is usually considered 
deleterious, and it is not generally recognized that ran- 
dom loss can suppress fluctuations, even below the Pois- 
son level. Thisis the atomic analog to intensity squeezing 
in optics 23|, 



24, 25 



26l ] . We show that three-body loss 
can be used to prepare small and well-defined numbers 
of atoms in each trap, ultimately enabling the study of 
collective excitations in mesoscopic ensembles. We trap 
a large number of dense mesoscopic ensembles in a lat- 
tice of microtraps which undergo rapid three-body decay. 
Through sensitive absorption imaging we measure the 
shot-to-shot distribution of atom numbers and find sub- 
Poissonian statistics for between 50 and 300 atoms per 
trap. The effects of residual imaging noise are greatly re- 
duced through the application of spatial correlation anal- 
ysis which exploits the lattice geometry and provides a 
way to isolate atom number fluctuations. Our results 



are in very good agreement with a model for stochastic 
three-body loss which takes into account the fluctuations. 

For ultracold gases in magnetic microtraps, inelastic 
density-dependent decay is the dominant loss process. 
In 87 Rb this is typically due to three-body recombina- 
tion 27,28, [29[ , whereby all three atoms are lost from 
the trap. As this depends on the probability of finding 
three atoms together, three-body recombination is a sen- 
sitive probe of density fluctuations and correlations in 
degenerate Bose-gases [30j, |3lj, |32J. This previous work 
involved the macroscopic evolution of the mean number 
of remaining atoms, which decays proportional to the 
mean square density. 

We are primarily interested in the fluctuations in the 
number of remaining atoms. We model this with the 
following master equation for the probability distribution 
P(N,t), 
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which is valid for any birth-death process with multiple 
reactions involving p bodies [33| . Here Nq is the initial 
mean atom number in a given trap, k p are the scaled rate 
constants and the step operator E p changes N — > N + p. 
Eq. (fTJ) is a set of coupled differential equations, one for 
each possible value of N. For small systems involving up 
to a few hundred atoms, these equations can be solved 
numerically to provide the full atom statistics (including 
fluctuations) as a function of time. 

In our experiments hi « and k% 3> fci. For a non- 
degenerate gas at temperature T in a harmonic trap, 
k 3 /N% = (2L 3 /V3)(mcj 2 /27rfc B T) 3 , where the mean trap 
frequency in our case is tD = 27r x 10.0±0.5 kHz. The 
three-body rate constant is L3 = 1.8(±0.5) x 10 -29 cm 6 /s 
for the T — = 2 hyperfine state of 87 Rb 31 1. 

For the mean and variance of the distribution we can 
obtain approximate analytic expressions. Following [33| 
we perform a system size expansion for Nq 3> 1 to obtain 
a linear Fokker-Planck equation and derive equations of 
motion for the moments. For combined one-body and 



2 




300 



C- 200 



100 





1 body 


3 body 




■ t * i * i i * 







0.05 0.1 

OD 



500 1000 1500 2000 2500 

t (ms) 



FIG. 1. (a) Subsection of an optical density image for 
t = 25 ms for a single realization of the experiment, (b) De- 
cay of the mean atom number in a selected site (highlighted 
in (a)). A fit to Eq. ((2| is shown (solid line) together with the 
corresponding one-body loss (dashed line) . The shaded region 
indicates the range of atom numbers in our data for all traps. 
The residuals from the fit and standard errors on the mea- 
surement of (N) are shown below, demonstrating agreement 
at the level of ±3 atoms. 



three-body loss the evolution of the mean fraction of re- 
maining atoms is 



V 



(N) _ cxp(-M) 

1% ~ y/1 + (k 3 /h) [1 - exp(-2M)j ' 



(2) 



We express the fluctuations in terms of the relative 
variance, or Fano factor, F = ({N 2 ) — (N) 2 )/{N), where 
the averages are taken over realizations (F = 1 for a 
Poisson distribution). The evolution of F in time can be 
written as a function of 77. This leads to the following 
differential equation: 



dF _ k 3 ri 2 (5F(ri) - 3) + ki{F(ri) - 1) 
dr\ r)(ki + k 3 rf) 



(3) 



In the case where three-body loss dominates we obtain 
the simple solution 



F(,) = | + , S ( Fo -|), 



(4) 



where Fq = F(j] = 1) is the initial Fano factor. As the 
atoms are lost from the trap the Fano factor asymptotes 
to a value of F — > 3/5, significantly below the Poissonian 
level F = 1. Correspondingly the memory of the initial 
Fano factor is lost very rapidly due to the fifth power of 77, 
in contrast to one-body loss where F = 1+t](Fq — 1). The 
result is easily generalized to an arbitrary p-body pro- 
cess yielding an asymptotic Fano factor F — > p/(2p — 1). 
The results of this simple analytic model are in excel- 
lent agreement with the numerical solution to Eq. |T]) for 
(N) > 10. 

Our experiment incorporates a two-dimensional lat- 
tice of optically resolvable magnetic microtraps produced 
by a magnetic film atom chip 34, [35|. We load a few 



thousand atoms into each of approximately 250 traps, 



and then evaporatively cool close to quantum degener- 
acy (temperature ~ 3 pK, phase-space-density ~ 0.3). 
After cooling, a few hundred atoms remain in each trap. 
Due to the small size of each trap the atomic density is 
high (» 2 x 10 14 cm- 3 ) and we observe rapid three-body 
loss, despite relatively few atoms per site. During the 
experiment we apply a fixed radio-frequency 'knife' (ef- 
fective trap depth ~ 35 pK) to ensure the temperature 
of each cloud does not vary. The knife counteracts any 
heating that may accompany the three-body loss and due 
to the high trap depth the role of heating-induced loss on 
the expected fluctuations is negligible. 

We image the in-situ distribution of atoms (Fig. QJl) 
using absorption imaging in reflection geometry with a 
circularly polarized probe laser aligned perpendicular to 
the chip surface [35]. The effective pixel size in the object 
plane is 3.2 pm, the optical resolution is 7.5 pm (Rayleigh 
criterion) and the lattice spacings are 22 pm and 36 pm. 
The exposure time is 0.15 ms and the saturation param- 
eter is s = 2 x 0.3 (double pass). In each run of the 
experiment we record an absorption image, a reference 
image taken without atoms and a stray light image, from 
which we compute an optical density image of the atomic 
distribution. Each image contains the center-most region 
of the loaded lattice and a surrounding background re- 
gion used to quantify the imaging noise. 

To measure the decay we hold the atoms for a variable 
time after evaporative cooling before taking the absorp- 
tion image. Our data are comprised of two sets. The 
first spans from t = ms to t = 880 ms with 40 intervals 
(selected on a power- law scale) and repeated 15 times 
(600 runs of the experiment). The second data set spans 
from t = 21 ms to t = 2.5 s, with 40 intervals and 19 
repeats (760 runs). From this we extract for each mi- 
crotrap (with index m) (i) the decay of the mean atom 
number (N m ) and (ii) the variance (N^) — (N m ) 2 . Each 
optical density image is aligned to the average to mini- 
mize the effect of jitter between shots. The atom number 
in each site and each image is found by a two-dimensional 
amplitude fit of a model shape function which minimizes 
the influence of imaging noise. The shape functions are 
obtained by Gaussian decomposition of the average op- 
tical density image, for each cell of the lattice. The ob- 
tained shapes are smooth peaked functions which repro- 
duce the observed absorption profiles (FigQJa)), account- 
ing for small distortions due to the underlying chip sur- 
face. Least-squares amplitude fitting is then performed 
on each image for 245 ensembles, with each fit including 
the 8 nearest neighbors to account for small overlapping 
areas. 

Figure [T^ shows a section of an optical density image 
for a hold time of 25 ms. The evolution of (N) for a 
selected trap is shown in Fig. Q}). A fit to the data with 
Eq. © yielding k 3 = 10.4 ± 0.4 s" 1 , ki = 0.52 ± 0.03 s" 1 
and Nq = 354 ± 4 is shown, together with the corre- 
sponding one-body decay 77 = exp(— kit). The shaded 
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FIG. 2. Atom number fluctuations measured for each of 245 
lattice sites during three-body decay. (a,b) Number distri- 
butions for one specific trap (m=138) at two hold times. 
Histograms correspond to 19 measurements and each bin is 
5 atoms wide. The lines indicate Gaussian fits to the data 
(solid), Poisson distributions (dashed), and combined Pois- 
son and detection noise contributions (dash-dotted), (c) The 
relative variance vs. {N m ) for each lattice site and for each 
hold time (points). Open circles indicate the measurements 
for the selected trap, with a fit (including detection noise) for 
a constant Fano factor F = 0.57 (solid line). Arrows highlight 
the two data points corresponding to the histograms (a,b). 



region indicates the full range of atom numbers in our 
data obtained by analyzing each trap individually. 

To quantify the fluctuations it is necessary to accu- 
rately calibrate the absorption cross-section. For this we 
compare for each trap individually the measured cloud 
temperature and three-body loss rate [35j |. to indepen- 
dently infer the atom number. In this way we de- 
termine an absorption cross section of (0.32 ± 0.05)cto 
((To = 3A 2 /2tt) which is in good agreement with the ex- 
pected cross-section of 0.31cto based on our imaging pa- 
rameters. The maximum optical depth for a trap con- 
taining 250 atoms is ~ 0.1. 

Figure [3] shows the measured atom number statistics 
for various hold times, corresponding to different mean 
atom numbers in each trap. A histogram of the fit- 
ted number of atoms in one specific trap for 19 rep- 
etitions of the experiment at t — 25 ms is shown in 
Fig. The measured (N) = 280 ± 3 and the variance 
is (iV 2 ) - (N) 2 = 140 ± 50, indicated by the Gaussian 
distribution (solid line). The distribution is significantly 
narrower than for a Poisson distribution (dashed line), 
providing a direct observation of sub-Poissonian num- 
ber statistics in our experiment. For longer hold times 
(Fig. f^b) the mean number of atoms decreases due to 
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FIG. 3. Lattice averaged Fano factor F as a function of the 
mean number of atoms (N). Horizontal lines correspond to 
F — 1 (dashed) and to F — 3/5 (dotted) for strong three- 
body loss. The solid line is a model including three-body and 
one-body loss terms. The shaded region indicates systematic 
uncertainties described in the text. The inset shows an ex- 
ample fluctuation correlation function (xW) f° r t = 25 ms. 



loss, however the observed distribution does not become 
significantly narrower. This is due to the added detec- 
tion noise contribution (dash-dotted line) which begins 
to dominate the observed fluctuations for (N) < 60. 

The same analysis is performed for each site and each 
hold time independently to obtain the site-resolved rela- 
tive variance as a function of the mean number of atoms. 
Fig. shows the results of 245 x 40 observations where 
each point is derived from 19 measurements. The ob- 
served fluctuations have two main contributions, atom 
noise with a constant F (Poisson noise is indicated by 
a dashed line) and a detection noise contribution corre- 
sponding to a fixed variance of 64 atoms 2 /trap/shot (dot- 
ted line). We find for N > 100 the vast majority of data 
points fall well below the combined variance for Poisson 
fluctuations (dash-dotted line) indicating F < 1. Inter- 
estingly the deviation from Poisson statistics is most ap- 
parent for small hold times (large (N) ) , indicating three- 
body loss also has a significant effect on the fluctuations 
before the end of the evaporative cooling stage. 

To account for detection noise and to investigate the 
sub-Poissonian noise over the full range of atom numbers 
in our experiment we perform spatial correlation analy- 
sis of the optical density images. Here wc benefit from 
the lattice geometry and separate various noise compo- 
nents based on their respective correlation length-scales 
to isolate the atom fluctuations in our data. 

We compute, for each optical density image the two- 
dimensional fluctuation correlation function Xi(^) = 
J(rii(x) — (rii(x)))(ni(x + 8) — (rii(x + 5)))d 2 x, which is 
then averaged over the realizations of the experiment 
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(indexed by i) for a given hold time (Fig. [3] inset). 
We model the observed spatial distribution by m(x) = 
c iYj m N i™P™( x )+ d i( x )' where N i<m &ndp m (x) arc the 
number of atoms and local shape function respectively for 
ensemble m, a ~ 1 accounts for correlated noise (due for 
example to probe frequency noise) and di(x) accounts 
for spatially uncorrclatcd imaging noise. The correlation 
function (x(<5)) shows several distinct features (Fig.[3]in- 
set). A narrow spike at 8 = (central red pixel) rep- 
resents the uncorrelated imaging noise. This sits on top 
of a broader peak (dark central feature) representing the 
fluctuations correlated over the length scale of approx- 
imately a single cloud which accounts for shot-to-shot 
fluctuations of the number of atoms within each trap. 
An array of neighboring peaks, spaced at the lattice pe- 
riod, represents the correlated noise across traps which 
we attribute to small fluctuations of the probe detuning. 

In the analysis of (x(<5)) f° r each hold time, we first 
subtract the calculated background-region correlation 
function and exclude the <5 = pixel spike to elimi- 
nate the uncorrelated imaging noise. We then fit two- 
dimensional Gaussian distributions to the central and 
neighboring correlation function peaks. The lattice- 
averaged Fano factor (weighted by (N m )) is given by 
F = (EmW ~ Y. m (N m ) 2 )/ E m (N m ). Neglecting the 
small overlap between neighboring shape functions and 
noting that the fluctuations of iVj, m are uncorrelated be- 
tween different traps, we obtain 



(5) 



where Xq and X\ are the fitted volumes of the central 
and neighboring peaks of (xW) respectively, Po is the fit- 
ted volume of the pre-averaged autocorrelation function 
peak J (iii(x))(ni(x+5))d 2 x, and (N) is the weighted av- 
erage atom number. Accounting for the overlap between 
neighboring shape functions yields a small correction fac- 
tor, which for our lattice geometry is < 1.1. 

Figure [3] shows the extracted Fano factor for two sepa- 
rately analyzed data sets as a function of (N) during the 
hold time. Horizontal lines correspond to the Poissonian 
limit F = 1 (dashed) and to the expected limit F = 3/5 
(dotted) for strong three-body decay. The data shows 
sub-Poissonian atom number fluctuations for (N) > 50 
up to 300 atoms per site. A fit over this range indi- 
cates a Fano factor of F = 0.53 with a standard devi- 
ation of ±0.08. We independently estimate a system- 
atic uncertainty of ±0.2 incorporating uncertainties in 
the absorption cross-section, background noise contribu- 
tion and the overlap between neighboring traps. The 
measured fluctuations are clearly below the Poissonian 
noise level (dashed-line) and are in good agreement with 
the theoretical expectation of F = 3/5 (dotted- line). For 
(N) < 50 one-body loss dominates and we expect F to 
increase to 1. The solid line is the result of Eq. ((3j) in- 
cluding both three-body and one-body loss terms. 



In conclusion, we have shown that normally- 
undesirable density dependent losses in small atomic en- 
sembles naturally lead to suppressed fluctuations of the 
absolute atom number to below Poissonian noise lev- 
els. By three-body decay it is possible to prepare hun- 
dreds of small and well-defined atomic ensembles con- 
sisting of tens to a few hundred atoms. We expect this 
to be an ideal system for the study of collective exci- 
tations produced for example via laser-excited Rydberg 
states for quantum information processing with neutral 
atoms [ljl 17, 18, Such ensembles also have de- 

sirable properties for generation of Schrodinger-cat-like 
states [19j , the study of spin-squeezing and as a resource 
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for quantum metrology using trapped atoms 
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